Antiferromagnetic Order in Pauli Limited Unconventional Superconductors 
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We develop a theory of the coexistence of superconductivity (SC) and antiferromagnetism (AFM) 
in CeCoIns. We show that in Pauli-limited nodal superconductors the nesting of the quasi-particle 
pockets induced by Zeeman pair-breaking leads to incommensurate AFM with the magnetic moment 
normal to the field. We compute the phase diagram and find a first order transition to the normal 
state at low temperatures, absence of normal state AFM, and coexistence of SC and AFM at high 
fields, in agreement with experiments. We also predict the existence of a new double-Q magnetic 
phase. 
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The interplay between antiferromagnetism and uncon- 
ventional superconductivity is one of the most intensely 
investigated topics in correlated electron systems. The 
115 family of heavy fermion compounds (CeMIns where 
M=Co,Rh, Ir) exhibits many salient features of this in- 
terplay. CeCoIn 5 is a clean layered singlet <i-wave su- 
perconductor whose upper critical field, H C 2, is deter- 
mined by the Zeeman splitting of the electronic states 
(Pauli limiting), rather than by the orbital motion of 
the Cooper pairs. This material has an unusual low- 
temperature phase at fields just below H c2 , separated 
by a second order transition line from the low field su- 
perconductor, and by a first order transition from the 
normal (N) state (H > H c -i) Initially this phase was 
conjectured [l[ to be the first realization of the Fuldc- 
Ferrell-Larkin-Ovchinnikov (FFLO) state, where Cooper 
pairs acquire a finite center of mass momentum to coun- 
teract the pair-breaking due to Pauli limiting, and the 
superconducting order parameter oscillates in real space. 

However, experiments present a more complex situa- 
tion. The high-field N-SC transition becomes first order 
at temperatures, T, higher than that of the onset of the 
new phase [l[. Specific heat and resistivity are anoma- 
lous just above H C 2, consistent with quantum critical be- 
havior due to antiferromagnetic fluctuations @, H[, yet 
AFM was not found upon suppression of SC by several 
means 0, D . Finally, a series of NMR and neutron 
scattering measurements with H parallel to the layers 
established the existence of static AFM in the high-field 
phase [7-12|. The ordered moment m _L H is modu- 
lated with an essentially field-independent wave-vector 
Q/(2tt) = (<7,g,0.5), (q ~ 0.44) jaEl. This unusual 
phase challenges our understanding of the connection be- 
tween SC and AFM. 

In this Letter we show that the near-perfect nesting 
of the quasiparticle pockets created by Zeeman pair- 
breaking near the nodal regions creates the conditions 
for the antiferromagnetic instability even when the nor- 
mal state Fermi surface is not nested. This AFM insta- 
bility is equivalent to equal spin pairing of Bogoliubov 



quasiparticles within each pocket. We compute the tran- 
sition temperature both analytically in the weak-coupling 
limit (WCL) and numerically for a two-dimensional (2D) 
model. The resulting coexistence of AFM and SC ap- 
pears in the low-temperature and high-field region of the 
phase diagram, while magnetic order does not appear in 
the normal state. Moreover, the N-SC transition becomes 
first order at a temperature above that where the mag- 
netic order first appears. These results are in agreement 
with the experimental observations on CeCoIns. The 
most favorable magnetic state in the WCL is a double-Q 
collinear structure at the incommensurate wave- vectors 
connecting opposite nodes, with m _L H, but a single-Q 
state may be stabilized in a part of the phase diagram. 

Previous theories focused on the field induced AFM 
and considered the single-Q state only. In some pro- 
posals the amplitude modulation of the superconducting 
order parameter, A(r), is essential and drives a spin den- 
sity wave (SDW) order that varies on the same length 
scale [lH, 13 • I n other theories the SDW instability is 
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not conditional on the oscillatory behavior of A(r) 
I3- In Refs. (lEIil 

the direction of the SDW magneti- 
zation was taken parallel to H ', and the phase diagram 
was obtained numerically for fixed (large) values of in- 
teractions in different channels. We find that the longi- 
tudinal antiferromagnetic susceptibility is gapped under 
Zeeman field, but the transverse susceptibility is diver- 
gent. Hence AFM with m _L H appears naturally. Lan- 
dau expansion of the free energy in a model with critical 
magnetic fluctuations also supports m _L H for nodal 
superconductors jT^j- It was recently suggested [IH that 
an enhancement of the density of states above the normal 
state value due to combined effect of vortices and strong 
Pauli paramagnetism leads to the SDW instability, with 
no FFLO-like modulation of A(r). In our theory it is the 
nesting between pockets, and therefore the joint density 
of states enhancement, that drives the magnetic transi- 
tion. Finally, we propose the existence of a double-Q spin 
structure that differs from previous theories and provides 
a test of our results. 
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FIG. 1: (Color online) Fermi surface and the nodal pockets 
generated by Zeeman field for the tight-binding model de- 
scribed in the text with Ao/t = 0.7, h/t = 0.3,0.7. feii and 
k± correspond to the directions parallel and perpendicular to 
the Fermi surface near each pocket. 



Since the main Fermi surface sheet of CeCoIns is quasi- 
cylindrical, we consider a 2D Hamiltonian for a d x 2_ y 2 
superconductor with a Zeeman term h = gfisH, 
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Here e k is the band dispersion and A k is the d x i_ y i or- 
der parameter. We select the spin quantization axis along 
H. A Bogoliubov transformation 7^ = u k c k ^ 



vie 



^[l + £fc/^fc],2| Ufc | 2 = 
A? brings the Hamilto- 



7*4 = UkCki+v* k c\ kx with 2\u k \ 
[l-e k /E k ], and E k = + 
nian to the diagonal form 

Ho = { l E * - h \ TkfTfct + \ E k + h] ll kl l-ki} ■ (1) 

k 

The first term describes the creation of spin-polarized 
pockets of Bogoliubov quasiparticles by magnetic field. 
These pockets are unstable to the formation of AFM, 
and we first elucidate this instability in the WCL. 

At low fields we linearize the dispersion in the near- 
nodal regions by taking e k ~ vpk± and A k « v&k\\ 
with k± (fcy) measured from the node normal (parallel) 
to the Fermi surface (see Fig. [1} and obtain elliptical 

Fermi pockets ^Jvpk^_ + v^k 2 = h. For a centrosym- 

metric system with pairs of nodes at momenta ±fe n j, the 
low field pockets are perfectly nested at Qi = 2k ni and 
Qi = —2k ni (i = 1,2, see Fig[T]), even though the Fermi 
surface in the normal state is not. This nesting drives 
the instabilities of the superconducting state. 

As 115 materials exhibit strong antiferromagnetic fluc- 
tuations, we assume that there is a residual interac- 
tion in this channel. We start by considering the sim- 
plest AFM with m _L H: a singlc-Q spiral m(r) = 



mq(cosQ • r,s'mQ • r,0), where Q £ ±Q» is one of 
the four wave- vectors that are perfectly nested. The mean 
field Hamiltonian for this spiral order is 



H A 



FM 



= -J 



-fc+QT + m Q c fc+Qt c *4 



(2) 



= J2k( c l+Qt c ki)/ N , J is the coupling con- 



where 

stant, (. . .) denotes a thermodynamic average, and N is 
the number of lattices sites. 

Since a Bogoliubov quasiparticle with a given spin at 
momentum k combines an electron state with the same 
spin at fe with a hole with opposite spin at — fe, the AFM 
interaction in the particle-hole channel is manifested as a 
spin triplet superconducting instability of the Bogoliubov 
quasiparticles within each pocket. Transforming Eq. ^fy 
into the basis of Bogoliubov states we find 



c fc+Qt Cfc J- 



-Uk+QVkll+Q^-M 



Note that we only need to consider the quasiparticles 
with the spin along the field (occupied states). The terms 
marked as irrelevant (irr.) include spin down Bogoliubov 
operators, and therefore do not have an expectation value 
since these excitations are gapped by h. Consider the 
vicinity of the node —Q/2. We have fe = —Q/2 + q, so 
that fc + Q = Q/2 + q, and we find 
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which describes the triplet pairing of "Bogoliubons" 
within the same pocket. Since u k+ qv k w sgn(fcy) in the 
linearized approximation for a d-wave gap, we obtain 



H = 



E -%+q 



+^™Q7s +(?t 7s_„ t s g n (9||) + H - c 



qi" 



(3) 



showing p-wave pairing symmetry. The sum in Eq.Q 
runs over wave-vectors q near the pocket center, i.e. 
\q\ <C Tr/a. This approximation is valid as long as Q 
is far from the commensurate wave vector (jt,tt). In that 
case, only the wave- vectors on the Fermi surface of the 
pocket centered at —Q/2 satisfy the nesting condition for 
the transferred momentum Q (sec FigJ5|). 
The critical temperature in the WCL is 



2it 2 vfva 
Jh 



(4) 



where uj c is the high energy cutoff for the magnetic 
interaction J. This result is valid for Tq <C h <C 
{w c ,Ao}, and describes the instability of the uniform 
SC towards the AFM modulation. It is likely that 
the exponentially small m at low H is destroyed by 
the orbital effect (see below). The order parameter 

has 



that appears at Tq, I]q s g n (9||)(7r 



Q/2+qt 7 Q/2-qt 
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a staggered triplet superconducting component, A 



--Q/2- 



-at -Q/2- 



Q 



?t ) with sgn((5* ) = sgn(g||), in ad- 



dition to the AFM rriQ |15|, [16|. This triplet compo- 
nent is an inevitable companion of Ao and ttiq since 
the Ginzburg-Landau expansion allows for the invari- 
ants of the form Aq • jtIqAq + Aq • tuqAq, linear in 
Aq, where Aq is the vector of three independent com- 
ponents of the triplet superconducting order parameter 
(2Aq = Aq + iAq). However, we expect Aq to be ex- 
traordinary fragile with respect to impurity scattering so 
it cannot be easily observed in CcCoIns. 

To go beyond WCL, we consider a mean field Hamil- 
tonian with both SC and AFM interactions, 



[ m Q c U c k+Qi + H.c] 



I An 



V 



J\m Q \ 2 , 
(5) 



with a tight-binding dispersion = 2t(cos k x + cos k y ) — 
\i. We fix n/t — 0.749 which sets the wave vectors con- 
necting the nodal points to Q = (±0.88tt, ±0.88tt). The 
order parameter A& = Aq.% with rjf, = cos k x — cos k y 
and Ao = VJ2k %( c -fe4- c fet)- Again we first consider a 
single- Q XY magnetic ordering. We diagonalize H and 
minimize the free energy to obtain the phase diagram as 
a function h and T for given interactions J and V. The 
normal state becomes unstable towards magnetic order- 
ing at J c = 3.6t. Consequently, we choose J = 3.5i 
to describe a system with strong AFM tendencies. The 
phase diagram for V = 3.0t in Fig. [3] shows the coex- 
istence of AFM and SC (Q-phase) below T Q < T* at 
high fields as well as the first order N to SC transition at 
< T c (h) < T* . The latter appears here due to purely 
Zeeman coupling to the field, dominant at high fields for 
Pauli- limited CeCoIns, but is known to persist in the 
presence of moderate orbital effects [19j . The phase di- 
agram is qualitatively the same for 0.9J C < J < J c and 
is robust against changes of V and the inclusion of the 
next-nearest neighbor hopping t'. This strong coupling 
result confirms the instability analytically found in the 
WCL [see Eq.flJ)], and shows that the phase diagram is 
ubiquitous for systems such as CeCoIns. 

The above analysis demonstrates the instability of a 
Pauli limited d-wave superconductor towards incommen- 
surate AFM. We now determine the most stable mag- 
netic order. If the nodes of the order parameter are 
located at (±l/y/2,±l/-/2)k F , there are four distinct 
wave-vectors, Qi = (±V2,±V2)kp that connect the 
nodal quasiparticle pockets. The most general magnetic 
structure in the SDW state is m(r) = m x (r) + im y (r) = 



m Q 3 



(recall the spin z-axis 



is along H). As CeCoIns has an easy c-axis (taken as spin 




FIG. 2: (Color online) Conditions for pairing of the Fermions 
at the nodal pockets due to the term c \ s+ Q ^Cki- a place- 
holder. Dotted line: Fermi surface, green dots: nodes. Left 
panel: If k is in the lower pocket, only the term with Qi = Q 
satisfies the nesting condition and leads to pairing. Right 
panel: If k is in the upper pocket, only Qi — —Q is relevant. 



vorable. Now consider each single-Q term c^Q^Cki- Re- 
call that we pair spin-up Bogoliubov quasiparticles and 
that the electron creation (annihilation) operator at fc 
is a linear combination of the Bogoliubov creation (an- 
nihilation) operator with the same momentum and spin, 
and an annihilation (creation) operator with the opposite 
spin and momentum. As Fig. [2] shows, for fc, for exam- 
ple, in the pocket centered around — fc n i, Eq.([2]) produces 
a logarithmically divergent pairing interaction only for 
Q = Qi = 2fc„i- In contrast, if fc is in the vicinity of fe n i, 
only Q = Qi = — 2fc„i is relevant for the AFM instabil- 
ity. Similar analysis applies to the other pair of pockets. 
Since Qi are far from (7r,7r), in the WCL the equations 
for the four Fourier components m,Q i decouple, look iden- 
tical, and hence |toqJ = mo for all Qi. The condition 
m y = immediately yields a doublc-Q collinear struc- 
ture, m(r) = m [cos (Qi ■ r + ip\) + cos (Q 2 ■ r + (p 2 )], 
with m\\c and simultaneous modulation of SDW along 
the [110] and [110] directions. Note that the stabiliza- 
tion of this structure is simply due to the counting of 
the minimal number of independent Fourier components 
of the order parameter that are needed to fully gap the 
Fermi pockets of Bogoliubov quasiparticles. The same 
counting leads to the single-Q collinear phase of Cr [20| . 

Beyond WCL, as states away from the Fermi sur- 
face participate in AFM, generally singlc-Q structures 
become advantageous 
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In CeCoIns elastic neutron 
scattering [9j finds a single-Q modulation with Q _L H 
for H along a nodal direction [110], and NMR line- 
shape at H > 10. 2T is consistent with a singlc-Q struc- 
ture [1, [l(J- At these high fields the large size of the 
Bogoliubov pockets leads to an intermediate or strong 
coupling regime. The first question is whether at lower 
fields, where pockets are small, our proposed double- Q 
structure is realized. NMR measurements of Ref. [10( 
show an anomalous broadening of the line shape for 
9.2T< H <10.2T and T = 70mK. Moreover, 9.2T is pre- 
cisely the field at which the amplitude of the ordered mo- 
ment of the single-Q structure extrapolates to zero, sug- 



x-axis) , the state with m y = is energetically more fa- gesting that another form of magnetic ordering persists in 
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FIG. 3: (Color online) Phase diagram obtained for V = 3t 
and J = 3.5t. Open and filled circles indicate first and second 
order phase transitions respectively. The red circles are the 
boundary between the SC and N states, while the blue circles 
are the boundary for the magnetically ordered phase. 



the 9.2T< H <10.2T region. The double-Q structure is 
a natural candidate for this intermediate phase because 
it produces a broad NMR line shape that is consistent 
with the experimental observation. Detailed neutron and 
NMR experiments will be able to test this hypothesis. 
Thermodynamic signatures of the transition from double 
to single-Q phase are expected to be weak. 

For if 1 1 nodes, the orbital coupling suppresses the 
SDW modulation with Q _L H at low fields. Within 
the semiclassical approximation the Dopplcr shift, 8E « 
vp ■ p s , where p s _L H is the momentum of the Cooper 
pairs, has opposite sign for the two pockets connected by 
Q, as vfi = —vf2- From Eq. (JlJ this expands one of the 
pockets and shrinks its counterpart, breaking the nesting 
condition. For some Fermi surface geometries an addi- 
tional field-dependent modulation may still allow SDW 
formation, similar to pnictides [22, 23 1 and exciton for- 
mation [24J. Generically however, AFM is favored when 



the typical Doppler shift 6E D ~ A ^H/H°^ b > fj, B H, 
where is the orbital critical field. In CeCoIns this 

gives fields H > Hp/H% b , where H P is the Pauli limit- 
ing field, above 20-40% of the observed H c2 (T = 0). The 
semiclassical approach is inadequate at high fields [HI 
where a more detailed analysis is required to determine 
whether the Q _L H state is stable. In weak coupling at 
low fields, SDW with Q\\H may exist since 8E vanishes 
for the quasiparticle pockets around the nodes aligned 
with the field [26| , but Tq is exponentially small. 

Anisotropic spin-spin interactions induced by the spin- 
orbit coupling may also favor ordering with Q _L H , but 
depend on specific models. An intriguing possibility is 
that this choice is due to putative FFLO state. In the 
presence of vortices the FFLO modulation is along the 
field at a wave vector determined by the field magnitude, 



and interferes destructively with the AFM at Q\\ H whose 
wave vector is fixed by the position of the SC nodes. 
Hence Q _L H is stabilized. Since generically the Tq(H) 
that we find differs from the onset of the FFLO mod- 
ulation, one expects two transitions from the supercon- 
ducting state, first to a double-Q AFM, followed by the 
FFLO+single-Q state. Further experiments will distin- 
guish between these scenarios. 

To summarize, our theoretical framework describes a 
generic instability of Pauli-limitcd unconventional super- 
conductors under Zeeman effect. This instability leads 
to incommensurate AFM at the wave- vectors connecting 
opposite nodal points. Our theory describes the existing 
experiments on CeCoIns and suggests studies to identify 
the double-Q collinear phase that we find to be favorable 
in the weak coupling regime, and which may exist in the 
intermediate field range. 
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